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Singular Phases of Seiberg-Witten Integrable Systems: 

Weak and Strong Coupling 

H.W.Braden 1 and A.Marshakov 2 

We consider the singular phases of the smooth finite-gap integrable systems arising in the context 
of Seiberg-Witten theory. These degenerate limits correspond to the weak and strong coupling 



p^ | regimes of SUSY gauge theories. The spectral curves in such limits acquire simpler forms: in 

most cases they become rational, and the corresponding expressions for coupling constants and 

qj | superpotentials can be computed explicitly. We verify that in accordance with the computations 

^"*^ ■ from quantum field theory, the weak-coupling limit gives rise to precisely the "trigonometric" family 

of Calogero-Moser and open Toda models, while the strong-coupling limit corresponds to the solitonic 

degenerations of the finite-gap solutions. The formulae arising provide some new insights into the 

fvq ' corresponding phenomena in SUSY gauge theories. Some open conjectures have been proven. 

>' 

o. 

^ ■ 1 Introduction 

Recent progress in understanding non-perturbative structures in supersymmctric (SUSY) gauge theories [f , 2] 
has also shed new insight upon the role of integrable structures in modern theoretical physics. Surprisingly, 
several integrable systems that were introduced and studied over the years as simplified models for quantum- 
mechanical and low-dimensional field theoretical problems have now reappeared as an important tool for un- 
derstanding physically interesting effects in almost realistic, effective gauge theories. Part of the reason for 
Oh. this (familiar from supergravity investigations) is that the scalars of say a D — 4 Af — 2 SUSY vector su- 

permultiplet are constrained to lie on a special Kahler manifold. With appropriate integrality constraints the 
cotangent bundle of such a manifold appears as the phase space of an algebraically completely integrable system 
(see, for example [3]). Yet an a priori argument for which integrable systems should arise, and what lies at 
the heart of this relationship between Seiberg-Witten (SW) theories [f, 2] and integrable systems [4] (see also 
[5, 6] and the books [7, 8] and references therein) still remains unanswered, and presents more questions than 
answers. Nonetheless this correspondence is accepted as providing a strong technical tool for the classification 
and computation of instanton effects in Af = 2 SUSY gauge theories [9, 10], with integrable systems yielding 
insights not yet understood in conventional gauge theoretic terms. Indeed one can straightforwardly present a 
dictionary for the correspondence between integrable systems and SW theories, but a majority of the "arrows" 
of the correspondence require further comment, not all of which is satisfactory. 

In what follows we shall study various singular limits of this correspondence. This will yield both a better 
physical insight and simpler expressions for the corresponding phenomena in the quantum gauge theories. The 
integrable systems arising in the context of SW theory are the well-known finite-gap solutions (see, for example 
[11]) associated to complex curves or Riemann surfaces of finite genus. Though comparatively simple complex 
manifolds, explicit formulae and constructions for these integrable systems and their solutions are often lacking. 
However, in certain singular limits, their structure simplifies drastically: this happens when one comes to the 
boundaries of the moduli space of their complex structures, or, in different terms, to the "boundary" values of 
vacua parameters of the gauge theory. These are the limits we shall study in this paper. The simplifications 
arising at these limits will enable us to derive explicit formulae and constructions, as well as provide certain 
insight into the corresponding phenomena in the quantum gauge theory. 

The boundaries of the moduli spaces of SW curves are characterised by either the vanishing or divergence of 
their period matrices Ty, which play the role of couplings in the effective gauge theories. The limits Ty — * oo 
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and Tij — ► correspond to weak and strong coupling phases of the gauge theory, and they are related by the 
modular transformation T <-> T~ x playing the role of S-duality [1]. Below we shall study the corresponding 
limits. Both limits appear rather naturally from the perspective of the integrable models. The weak-coupling 
limit of the gauge theory, when instanton contributions are exponentially suppressed and can be neglected, is 
almost trivial from the field theory perspective. The corresponding limit on the integrable systems side is a 
"simplification" of the interparticle potential: the non-periodic (or open) Toda chain is a limit of the periodic 
Toda "molecule" ; the trigonometric Calogero-Moser-Sutherland model is a degeneration of the elliptic Calogero- 
Moser model, and so on. In these cases (and now in fact only in these cases) can the SW integral formulae be 
practically derived from perturbative computations in the Af = 2 gauge theories. 

The strong-coupling limit of SW theories corresponds to another limit of the integrable models, that of 
smooth finite-gap solutions degenerating into solitons [12]. Physically this is the weakly-coupled limit of a dual 
theory in which one expects monopole condensation, breaking symmetry down to Af = 1, and confinement 
[1]. As we shall see this is in precise correspondence with the properties of the solitonic solutions of the SW 
integrable systems, which correspond to particular values of the integrals of motions. The corresponding vacuum 
expectation values (VEV's) are restricted to certain points in the moduli space, which can be treated as extrema 
of an Af = 1 supcrpotcntial. 

The paper is organised as follows. In section 2 we begin with a discussion of the weak-coupling regime 
in Af — 2 SUSY gauge theories, presenting some intuitive motivation for why this can be directly related to 
the integrable systems of the Toda or Calogero-Moser families. Here we also review general properties of the 
corresponding finite-gap solutions. In section 3 we discuss in detail the perturbative limit, giving rise to the open 
Toda chain family. We present explicit formulas for the generating functions for the open Toda chain (Toda 
molecule) obtained along these line and discuss their relations with matrix models and duality in integrable 
systems. We discuss in a similar way the perturbative limit of elliptic Calogero-Moser models. In section 4 we 
turn to the solitonic or strong-coupling limit of the finite-gap integrable systems and derive the explicit form 
of solution in terms of the Baker-Akhiezer functions. This is the most physically attractive phase of the SW 
theories, corresponding to the confinement in Af = 1 gauge theories, and we show that the phase of the solitons 
is related to the string tensions in the confining phase. We also prove (in Appendix B) a conjecture of Edelstein 
and Mas related to this strong coupling regime. We conclude with a discussion. 

2 J\f = 2 SUSY gauge theories and integrable systems 

We begin in this section by reviewing various aspects of the correspondence between Af = 2 SUSY gauge theories 
and integrable systems that we shall need in our later discussion. Throughout we will limit our attention to 
those models arising in the context of either D = 4 pure Af = 2 SU(N) gauge theories, or those with adjoint 
matter. These are listed in Figure 1. Extensions to other gauge groups and matter multiplcts are possible, 
though will not be our focus here. Indeed, for concreteness we will describe the N = 2 cases of these models in 
some detail. 

d = 5 Ruijsenaars-Schneider — ► 'Relativistic' Toda 

I I 

d = 4 Elliptic Calogero-Moser — ► Periodic Toda 

Af = 2 SYM with adjoint matter Pure Af =2 SYM 

I Perturbative Limits J. 

Calogero-Moser-Sutherland Open Toda 

Figure 1. 

We will begin the section with some intuitive motivation as to why the weak-coupling regime in Af = 2 SUSY 
gauge theories can be directly related to the integrable systems of the Toda or Calogero-Moser families. Then 
we review general properties of the corresponding finite-gap solutions and discuss in detail their perturbative 
limit, giving rise to the open Toda chain family. 

Although much of what we say in this section is by way of review we will also present several new technical 
details concerning the relation between the Calogero-Moser and Ruijsenaars-Schneider spectral curves. The 



books [7, 8] contain both recent reviews and general references for aspects of this correspondence not touched 
upon here. 

2.1 Motivations: Perturbative Calculations and N = 2 Examples 

The connections between SW theories and integrable systems can be discussed even at the perturbative level, 
where Af = 2 SUSY effective actions are completely determined by the 1-loop contributions. In the most well- 
known example, that of pure N = 2 gauge theory with SU(N) gauge group, the scalar field <fr = \\<&ij\\ of the 
J\f = 2 vector supcrmultiplet may acquire nonzero VEV's «I> = diag(ai, . . . , a^r) in cxtremising the potential 
Tr[<I>, 3r] 2 . This (generically) breaks the SU(N) gauge group down to the maximal compact [/(l)^ -1 . The 
masses of VF-bosons and their superpartners are (classically) given by m^ = a%j = a» — aj because of the 

interaction term ([A^, *&]ij) = \A l Hai — aj)) ■ By a simple technical trick these masses can be uniformly 
expressed in terms of the generating polynomial 



w = P N (X) = det (A *) = H(\ - ai ), (1) 

via the contour integral 

I ds ,pcrt = I XdlogW^J) \d\ogP N (X). (2) 






For a particular "figure-of-eight" like contour CV, around the roots A — a% and A = aj this may be computed 
via the residue formula. This means that the contour integrals in the complex A-planc (2) for such contours 
around these singular points (the roots of the polynomial (1)) give one set of the BPS masses in the SW theory. 
Another set of masses (monopoles) may be associated to dual contours starting and ending at the points ai 
where dS pcrt has singularities. Although such integrals are divergent, the result of such an integration says 
that the monopole masses are proportional to the masses m^ , multiplied by the inverse square of the coupling. 
The divergences are absorbed by a renormalisation of the coupling. The integrals over the contours on the 
marked plane we are dealing with here are, as we shall later see, to be viewed in terms of contour integrals on 
a degenerate Riemann surface. 

Now in M = 2 perturbation theory the effective action (or prepotcntial) T and the set of effective couplings 
Tjk, related by —iirTjk = d .* , are determined by the 1-loop diagrams. This results in the logarithmic term 

T »r = — v iog (mass)2 = ^iog4i ,. 

^ fe 2tti ^ g A 2 2m g A 2 ' ( 3 ) 

(massesM)jfc 

where the scale parameter A = Aqcd so introduced may be related to the bare coupling r. In the perturbative 
weak-coupling limit of the SW construction this is all one has. Instanton contributions to the prepotcntial which 
are proportional to powers of A 2N are suppressed, and one keeps only the terms proportional to log A or the 
bare coupling r. 

It remains to connect this discussion with integrable systems. The connection comes by interpreting the 
differential <i5' pert with the generating differential of a Hamiltonian system. Remarkably, this interpretation 
extends beyond the perturbative regime! Lets consider the simplest SW theory [1], the SU(2) pure gauge 
theory, where eq. (1) turns into 

w = A 2 - h. (4) 

Here h = iTr<I> 2 = ^a 2 and the masses (2) are now defined by the contour integrals of 

, ,, XdX XdX 

dS - Xdl0gw ^x-Th + x-W = h (5) 

Now, cqs. (4) and (5) can be interpreted as an integration of a simple dynamical system, the SL{2) open Toda 
chain (or Liouville model) with Hamiltonian 

H= l -{p i +e^). (6) 

This has solution e q = aj cosh(a[t — to]), where a = y2H. One verifies that 

p 2 dp apda , . 

-o 2 = P d Q + ~2 2 (') 

p z — a A p A — a A 



and so upon integrating over a trajectory § \_ v -± — § pdq — 2iria. But with the co-ordinate w = e 2q , 
momentum p = X and Hamiltonian (energy) H — h this is essentially (5). That is, integration of the canonical 
differential dS = 2pdq over the trajectories of the Toda chain solutions gives rise, for various (complexified) 
trajectories, to the BPS masses in the SW theory. 

This is actually a general rule. The perturbative mass spectrum and effective couplings of the Af = 2 theories 
of the "SW family" appear as actions when integrating the canonical differential of "open" or trigonometric 
families of integrable systems, the open Toda chain, the trigonometric Calogero-Moser and the trigonometric 
Ruijsenaars-Schneidcr systems. For pure gauge theories we have the iV-particle (open) Toda chain with (ratio- 
nal) curves given by eq. (1) and generating differential given by (2). The inclusion of adjoint matter is associated 
to the trigonometric Calogero-Moser-Sutherland model where 

W = ^ ){X) dS = X^. (8) 

Five dimensional gauge theories whose D — 4 reductions have adjoint matter correspond to the trigonometric 
Ruijsenaars-Schneider system, where now 



P ( N RS) (Xe- 2 ^) 



dS = \ogX — . (9) 

w 



In each of these cases Pjy and P N are appropriate polynomials that we shall further describe in due course. 
It is easy to see that (the perturbative) spectra are given by the general formula [13] 

M = aij® — ® — jr— , neZ. (10) 

H K 

In addition to the Higgs part a^ this contains the Kaluza-Klcin (KK) modes ^ and the KK modes for fields 
with "e-shifted" boundary conditions. The effective couplings arc defined by almost the same formula as (3) 



T-si £ ■*£- (id 

(masses M)jj 

i.e. by the sum of all the logarithms of spectrum (10) giving contribution for a particular 1-loop diagram for 
T 

A priori nothing so simple can be said about the spectrum and structure of the theory at strong couplings. 
However the structure of the "strongly-coupled" phase is, as we will see below using the correspondence with 
integrable systems, similar in many respects to the weak-coupling regime, and may be described by explicitly 
computable expressions. 

2.2 Generalities: finite-gap or Hitchin systems 

Let us now consider the general setting. Our analysis so far has been built upon a (polynomial or rational) 
relation R(X, w) = 0, a differential dS and BPS data given by integrating dS over various contours. We will now 
consider i?(A, w) — as the equation describing the spectral curve of an integrable system. The differential, 
contours and symplectic form may be canonically described in terms of this integrable system. At heart will be 
the choice of integrable system. For the SW theories under discussion these are given by finite gap integrable 
systems of particles of a particular kind, the so-called Hitchin systems [14, 15]. We shall now describe some of 
these general features and give explicit descriptions of the Hitchin systems of relevance to us. The remarkable 
feature of the SW integrable systems correspondence is that non-perturbative aspects of the field theory are 
incorporated into this construction [5, 6]. 

Replacing the relation R(X,w) = of (1), (8) and (9) one now has the Lax equation of a spectral curve, 

dct(A - C{z)) = 0. (12) 

Here the Lax operator C(z) is a matrix, defined on some base curve £o on which the spectral parameter z 
lies. For us this base curve is usually a torus (for elliptic models) or a sphere with punctures (for rational or 
trigonometric 3 models). The two further ingredients were the generating differential dS and the contours Cij. 

3 In perturbative examples it is parameterized by w = e z . 



The generating differential can now be denned by 

dS = Xdz 

( 13 ) 
^moduii^ = holomorphic differential. 

The action variables are given by the Seiberg-Witten contour integrals over half of the independent contours 

dS, (14) 

'A 

or 

a D = I dS. (15) 

where A = {Ai, . . . , A g } and B = {B\, . . . , B g } is a standard homology basis of H 1 (S). 
The holomorphic variation in (13) is crucial to the definition of the symplectic form 

= 5dS\j = da A dz(j) = dp A dq. (16) 

The variation here is to be computed at the divisor 7 = {71, . . . , 7^} of the poles of the Baker- Akhiezer function 
(which is determined by C{z)) which are co-ordinates on E s . Let {dcoi} be the the set of canonically normalised 
holomorphic differentials §. dujj — 5ij. The variation of dS is understood as a total external derivative on 
Moduli,, kS», 

5dS = 5(Xdz) = <5 mo duiiA A dz + dX A dz. (17) 

For this one must choose a connection on the bundle over moduli space [16] such that 5 mo duiiZ = 0, i.e. the 
parameter on the base curve z is covariantly constant. In practice this means that in the equation on two 
variables A and z which defines the spectral curve E, we consider z as an independent variable, with the 
variable A depending on moduli through eq. (12). Since on a one-dimensional curve £ for any differentials one 
has dX A dz = 0, from (17) we finally get 

dX ddS 

SdS = <5moduiiA A dz = Scii A — — dz — Sdi A — — = Scii A dcui. (18) 

OCli Odi 

Finally we introduce co-ordinates z on the Jacobian by 

9 ,-Y, 



p 


S pj 


1 dui + zVt 


-/] ^ + z ; 


7o 


j=l J 7o 



(19) 



with z\ = z\ (70 ) . Then 



ft = 2_, SdS(jk) — Sa-i A 2_. duiiilk) — 8a,i A 5zi- (20) 



(19) 

fe=i fe=i 



Our discussion thus far is general, predicated only upon a Lax operator and its attendant Baker- Akhiezer 
function. We must now describe the Lax operators for the class of models under consideration. These will 
be particular examples of Hitchin systems. For these the Lax operator in (12) C(z) = $ is considered as a 
meromorphic matrix-valued function (or better, a l-diffcrcntial) on the base curve S [11, 14, 15] satisfying 



<9$ + L4,$] = ^J (a) (5^(P-P Q ). 



(21) 



The right hand side serve as sources, and j( a ' are matrices whose structure (given shortly for a single puncture) 
is particularly simple for the su(n) theories. Thus $ is holomorphic in the complex structure determined by A 
on the punctured curve T,o/{P a }. The invariants of A can be thought of as co-ordinates (one commuting set 
of variables) while the invariants of $ as hamiltonians (another commuting set of variables) of an intcgrablc 
system. The most general features of these "holomorphic" finite-gap [11] or Hitchin [14, 15] systems are: 

• The spectral curve S g covers some base curve So (typically of genus go = 0, 1) and the g moduli of 
the cover are viewed as distinguished with the moduli of the base curve viewed as "fixed" or "external" 
parameters" . This set of distinguished moduli arc to be the constants of motion of our integrable system. 



• In the general context of finite-gap integrable systems one may view this set-up in a different way. A 
generic complex curve depends on 'ig — 3 parameters describing the complex structure. The integrable 
system is described [11, 16] by two meromorphic differentials dX and dz on the curve. By adding holo- 
morphic differentials these may be taken as having vanishing A periods. Fixing the B periods specifies 
the integrable system. Loosely, dX and dz are only defined up to multiples, and fixing the B periods 
of dX gives g — 1 constraints (taking account of the scaling freedom), while the B periods of dz give 
a further g — 2 constraints (now allowing dz — ► adz + (3dX). Altogether we come to a system with 
(3g — 3) — (g — 1) — (g — 2) = g parameters, the genus of S. 

• This construction (when g a = 0, 1) implies certain linear relations on the period matrix T^: 

Yl T v =T = fixed - (22) 

3 

This comes from the possibility to choose the homology basis of £ such that the A and B cycles "project" 
to the A and B cycles on the base torus So- (A rational base curve can be thought of as degenerated 
torus.) Then 



> Oj = y <b Xdz = <p ( y XA dz = y Cu &> dz — h\ i> dz = hi 



(23) 



and 

2_] af = 2_. f ^dz = (f> ( yj Xi 1 dz = \_] At f dz = hi (f> dz = hiT. (24) 



Thus, for any j, 



giving (22). 



\r- T ST 9 ®? d ST d (24) dhi (23) 

r i J J j J 



These properties we will be crucial in our discussion of the tau-functions (basically the Ricmann theta-functions 
on the corresponding spectral curves E) for these finite-gap integrable systems. 

Let us now give three particular examples of the Hitchin integrable models arising in the context of SW 
theory: the Toda chain, the elliptic Calogero-Moser model and the Ruijsennars-Schneider model. These are 
integrable systems we encountered in the perturbative discussion earlier. 

Toda chain. The Toda chain system [17] is a system of N particles with nearest neighbour exponential 

interactions: 

^l= Pi ^^eft+i-w.eft-fli-i. (26) 

dt dt 

It is an integrable system, with N Poisson-commuting Hamiltonians, h\ = ^2pi = P, hi = ^2 (^pf + e 9i_9i_1 ) = 
E, etc. We may have either an open chain (qa = —00, qN+i = 00) or a periodic system (of "period" N: qi+N = qi 
and p i+ 7v = Pi)- The periodic problem can be derived by reduction of the infinite-dimensional system of particles 
on the line by aid of two commuting operators: the Lax operator C (for the auxiliary linear problem for (26)) 

A*„ = J2 ZnkVk = e^"+i-9")*„ +1 +p n ^ n + e^"-'"- 1 ^-! (27) 

k 

and a second operator, the monodromy or shift operator in the discrete variable (the particle number) 

Tq n = q n +N Tp n = p n+N T*„ = ^ n +N- (28) 

The existence of a common spectrum for these two operators 

£* = A* T* = w* [£, T] = (29) 

means that there is a relation between them V(C,T) = 0. This in fact is the equation of the spectral curve S: 
V(X,w)=0. 



One way to get explicit form of the spectral curve equation is to rewrite the Lax operator (27) in the basis 
of the T-operator eigenfunctions. 4 This then becomes a N x N matrix, 



£ = £ JO H =p-H , ^ 

simple OL 



V" e a ^(E a + E_ a )+w- 1 e- a ^E_ 






c;a 



Pi 
5(92-91) 



C:2 



(92-91) 





.. we^ qi ~ 


-Qn) 


\ 


Vl 


22(93-92) 









(93-92) 


P3 










V 



1 1 

— e2 



(gi-9w) 











PN 



(30) 



/ 



Here the sum is over the simple roots, and olq = — ^simple a a IS m i n us the highest root. The construction 
depends explicitly on the eigenvalue w of the shift operator (28). This is the spectral parameter which is defined 
on the cylinder So = CP 1 \ {0,oo}. The eigenvalues of the Lax operator (30) are defined from the spectral 
equation 

r(\, w )=dct N (x-c TC ( w ))=o. (3i) 

Substituting the explicit expression (30) into (31), one obtains 



iV 



which for N — 2 is 



V(\,w)=P N (\)-w = 0, 

w 



1 



P N (\) = \ N + J2(-V k h k \ N - k , 



fe=i 



w+- = \ z - ^ +p 2 )X +p lP2 - (e q2 - qi + e qi - q2 ). 



The generating 1-form dS = A— indeed satisfies (13) 

SmoduhdS = 5 moduii dS\ w _ const = (<5 mo duiiA) — = ... = V" Sh k = holomorphic. 

w P N (\) w ^— ' y 



where 



y 2 = 



( w ~^) =^(A)-4, P' N (\) = 



dP, 



N 



OX 



(32) 
(33) 

(34) 
(35) 



h^— const 



By the gauge transformation [7y = u l £jj with w = v the Lax operator (30) can be brought to another familiar 
form 



C TC {w) - £ TC » = CT- 1 r rcr (t£;)l7 



p H + t," 1 e-«"-9S_ ao + £ simple a e a ^a + « e" ^*,, + £ 



/ Pi 

^52(92-91) 





apl 
I e i(Q2-9i) Q 

Ig2(93-92) 



./simple CK 



r.Oi 



«£-, 



w 2 



P2 
2(93-92) 



P:i 



V t 



I e 2(gi-9iv) 



^gKgi-gw) \ 






PA' 



(36) 



/ 





In this form [18] it clearly satisfies the 9-equation (21) on a cylinder with trivial gauge connection 

d „ 



B v C™{v) = ££tc {v) 



e-to-QE. 



-a, + E e0t ' 9E « ] *( P «>) 
simple Oi 



(37) 



*°-<*E ao + J2 e aq E- a \5(Poo). 

simple OL 



4 If we had chosen to work instead with C, which is a second-order difference operator, we would come to the Faddccv-Takhtajan 
2x2 formalism of Toda chains. 



The Toda chain coupling A (equal to the mass scale Aqcd in the SW theory) can be restored by rcscaling 
A — ► A/ A, w — * w/A N and the hamiltonians hk — > j*k- Then equation (32) acquires the form of 

K2N N 

w + =P A r(A) = A Ar + y(-l) fe / lfe A Ar - fe (38) 

fc=l 

proposed in the context of SW theory [2] . 

Elliptic Calogero-Moser model. The N x N matrix Lax operator (36) (and, therefore, (30)) can be 
thought of as a "degenerate" case of the Lax operator for the TV-particle Calogero-Moser system [19] 

C CM (z) =p-H + £ a F(q-a|«)£ a 

I p! F(q 1 -q 2 \z) ... F(q 1 -q N \z)\ 

F(q 2 - qi \z) V 2 ... F(q 2 -q N \z) (39) 

\F(q N -q 1 \z) F(q N -q 2 \z) ... p N J 

The sum here is now over all of the roots and the base curve So is a torus instead of a cylinder. The matrix 

elements 

, , .. o{q + z) , 

F ^ Z) = m ^j (40) 

are dchncd in terms of the Weierstrass sigma-functions a(z). Equivalently a(z) = 2ioe D ^~ g ,, z 2 where 0*(z) = 

9[ 1 ](z) = 9\{z) is the (only) odd Jacobi theta-function. The modulus of the elliptic curve So is r and we have 
the marked point z = at which x = p(z) = oo, y = ^p'(z) — oo. The canonical holomorphic 1-differential 
on So is dz = j?.. The Lax operator (39) corresponds to a completely integrable system with Hamiltonians 

h l = J2iPi> h 2 = J2i<j(PiP3 ~ m2 p(qi - Qj)), h 3 = J2 t<1< k(PiPjPk + ■■■), etC. 

-•CM 



From (12), (39) it follows that the spectral curve Yr for the A^-particle Calogero-Moser system 

dct(A-£^(z))=0 (41) 



covers N times the base elliptic curve. These spectral curves are very special: in general the genus of a curve 
dchncd by an N x N matrix grows as N 2 while the curve (41) has genus g = N. For N — 2 we have 

= A 2 - (p! + p 2 )\ + p lP2 - m 2 (p(q 2 -qi)- p(z)) . (42) 

The BPS masses a and & D are now related to the periods of generating 1-diffcrcntial 

dS CM = 2\dz = A— (43) 

y 

along the non-contractable contours on T, CM . 

In order to recover the Toda-chain system, one has to take the double-scaling limit [20], when m and — it 
both go to infinity and 

% - Qj -> ^(* ~ - ? ') logm + (* ~ q ^ ^ 

so that the dimensionlcss coupling r gets substituted by a dimensionful parameter A^ ~ m N e t7!T . The idea 
is to separate the pairwise interacting particles far away from each other and to adjust the coupling constant 
simultaneously in such a way, that only the interaction of nearest particles survives (and turns in an exponential) . 
This limit is described in more detail in Appendix A. In this limit, the elliptic curve degenerates into a cylinder 
with coordinate w — e z e nT so that dz — > — and 

w 

dS CM _ dS TC = A—. (45) 

w 

The Lax operator of the Calogero system turns into that of the A^-periodic Toda chain (30): 

C CM {z)dz^C TC (w)— (46) 

w 



m 



,2 V^ 



and the spectral curve acquires the form (31). 

One further remark is in order concerning the dependence of the elliptic Calogero-Moser model on the 
coupling constant, or cquivalently the mass of the adjoint hypermultiplct m. The exact equivalence between 
the Calogero-Moser and KP theories is usually considered only when m 2 = 2 (or at least m 2 = n(n + 1) for 
integer n). This restriction is however only essential when we consider the "Lax" equation for the "first" KP 
time, the x-variable, 

&_ 

dx 

highest pole). However, for the Calogero-Moser equations themselves (^pr = m 2 Si<j p'ilij) an d the "higher" 
hamiltonian equations) there is no such restriction on the coupling constant. One can always set m 2 — 2 by a 
simple rescaling of the time-variables t — > wit. The point is that in achieving the canonical forms of the KP and 
KdV equations all of the possible scalings have already been employed: consequently the coupling constant is 
restricted. 

Calogero-Moser spectral curve from the elliptic Ruijsenaars. The easiest and most general way to 
look at Seiberg-Witten theories with adjoint matter is to consider the elliptic Ruijsenaars-Schneider model [13]. 
The Lax operator of the elliptic Ruijsenaars model has the form [21] 



ih . 2 . - ^p{x- qi )U = E^. (47) 

This has solution with simple poles \& ~ J^ — ^— + . . . only when m 2 = 2 (the condition of cancellation of the 

'£&. -^2- 



k=£i 



In the trigonometric limit this turns into 



ePi sinh( gij +, + e)sinh(e) ^^tt L.^ ( } 

13 smh(^+e)smh(z + e)' 11 y sinh 2 (q lk ) 

Introducing Vi — e 2qi , ( — e 2z and q = e 2e one finds that 

qvi - Vj qC,-l C^°o qvi — Vj \( J £f Vi — Vj y ' 

Often only the leading term in (50) is taken as an expression for the Lax operator of the trigonometric Ruijsenaars 
system. 

Using (48) the spectral curve equation for the elliptic Ruijsenaars can be written as 5 

JV 

, T , . * — * ~w—w ~ rr I fii -• A- 7. A- f 1 rr I (= 1 

(51) 



det (A - C R (z)) = Y \ N - k (-) k Y X\e p det g fc + z + e W e ] 



fe=o \i\=kiei 



i,j£l 



where, using Wick's theorem for the correlation functions of free fermions [22, 21, 13], for the determinants one 
gets 

, a(q ij+ z+e)a(e) _ , 6,(q t] +z+e)6,(e) _ , e^qj-qj+z) 6,(z)6,(e) 

UeLfcxfe a(q i3 +e)a(z+e) - aCZ kxk 0„ ( 9 , 3+e )0„ ( z+e ) - " aCl kxk e,(q t - qj )8,(z) 6,(z+e) 



qj=q,-e 

k _ ,,,,,, a li 



^fff) det fcxfe (v(^fe)), = (%SSS i ) <n ie /«iWfe)>z 



k 

n.<j"" 

W t , J e,(q l -q ] ) ' 6,(z) 



( e,(z)e,(e) y rii<je.(gi-gj)rii<je.(gi-gj) e,(E <?»-£ <?»+z) 

V M*+<0 ) 

k 



(52) 



( e,(z) Y IIi<jg.(gij) e,(z+ke) 

\e,(z+e)) 



(z+£ 
,{z) 



Using (52) and introducing A = A " g ,) , one finally gets for (51) 



Rf~\\ - i 9 *( z ) 



det (\-C R (z)) = . 
nxn k y " \$*(z + e) 



N 



_k=0 

where the hamiltonians hu are 



\^X N ~ k h *( z + ke ) 

^ k 8*(z) 



(53) 



(_\ku, = V TT r P ' n^jg/fl.tg.j) _ V T-r p t T-r 6,(q tj ) 2 

{ ) ri k - 2^\i\=klUei e l[^ jeI 8,( qij +e) ~ 2^\i\=klUei e \.U<jei e,(q ij +e)e t (-q ij +e) ( 54 ) 
= z2\i\=kU. ie i e ' U. i<3 ei KW p(0-p(9ij) 



5 This is a particular case of generic determinant formulas (147), considered in more detail below. 



where we used 

8*(q + e)d*(q - e) u(g + e)a(q - e) 



= p(e) - P(<?)- 



0*(g) 2 0*(e) 2 ^(?)V(e) 2 

For JV = 2 we obtain 

det 2x2 (A -£*(.)) =A 2 -A( e ^+e^) +e ^^^^g ) 

= A 2 - Acr(e) (e»* + e P2 ) Vp(e) - p(gi 2 ) + e^+^a 2 (e) (p(e) - p(z + e)) . 

The spectral curve of the elliptic Calogero-Moser model (39) arises in the e — ► limit of (53). Indeed, in the 
limit e = mR, pi — ► Rpi and R — > 0, one recovers from (48) 

££ - e p * - e*« [] ( T(mi?)Vp(mi?)-pfe fc )- 1 + R Pi + 0(R 2 ), 



(56) 



C R = e* n ««^1 l^t'tw^i - * ' ™^4 + 0(tf) i * ,. 

^f (J{qi 3 + mR)<j{z + mR) a(q l j)(T{z) 



Thus 



43 - % + « fan + ro(l - %) ^ } ^ ) + 0(#) = % + AC« M + 0(^ 2 ) (57) 



(58) 



and one obtains 

dctjvxAr (A - C R (z)) = detjvxiv (A - 1 - R£ CM (z) + 0(R 2 )) 

= R N dct NxN (^-C CM (z))+0(R N + 1 ). 

Provided the Ruijsenaars parameter scales as ^r- — > A CM then in the limit R — ► we obtain exactly the 
spectral equation (41). 

Lets consider these for the N — 2 example. In this case one has from (55) 

det 2x2 (A - £*(*)) =(A-l) 2 -mi?(A-l)( Pl +p 2 ) 

+ (mR) 2 (Ap(g 12 ) - |A(p 2 + p 2 ) + ±( Pl + p 2 ) 2 - p(*)) + 0{R?) 

After the substitutions A — 1 — > i?A', one finally gets from (59) 

dot (i?A' + 1 - L R {z)) = R 2 (A' 2 - A'(pi + pa) + P1P2 + ™ 2 p(<Zi 2 ) - m 2 p(z)) + G(i? 3 ) ( 60 ) 

which is the N — 2 elliptic Calogero-Moser curve up to R N+1 = R 3 terms. On the other hand the N = 2 
spectral curve (59) can be rewritten in the form (53) 

det (A - £*(*)) = ( °W \ 2 (x 2 A g(z + C) /if + ^e^^ (61) 

2x2 v /y Vcr(z + e)y V cr(z) fj(^) / 

where 

X = \^1±J1, h R = (e^+encT(eWp(e)-p(q 12 ). (62) 

<t(z) 

With pi — > epi then in the limit e = mR — > this may be expressed as 



A 2 - \^±f h R + £^+H£i e Pi+P2 = (X _ 1)2 _ e (A - i)(^fM + 2 £^) 

+ e 2 (A/if M + (2 - A)(^f - ^f M ) + |(A - l)(ftf M ) 2 ) + G( £ 3 ). 



(63) 



Here we have used that h R = 2 + e/if M + e 2 (5 (/if M ) 2 - /if M ) . Finally after the substitution A- 1 -*■ eA', one 
arrives at 

A 2 - A^i /if + £^+ 2 i) ePl + P2 = e 2 ^,2 _ Va cm + h CM + {h CM _ 2A')^M + ^gi) + G(e 3 ) 

= e 2 *(«"&) (A' 2 - A'/if M + /if M ) + G(e 3 ). 
One can also rewrite (64) as 

A' 2 - A'/if M + /if M + (/if M - 2A') 2-^) + <£gl 



A'-^>) -( A '-^)/if--^ + ^-/if- (65) 

A 2 - Aftf a/ + frf W - p(z) 
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which is a common representation for the N — 2 Calogero curve (see, for example formula (108) below, with 
A = A' - £&, x = p( z ) and h$ M = -h). 



Returning to the general setting, instead of (64) one gets 



N 



\N-k h 0*(z + ke) 



"£~x N - k h k 



JV' 



{■- >>y>.J2\' N - k h% M + o(e N+1 ) 



<?.(*) 



(66) 



fe=0 ■ ' v ' k=0 

which is the D'Hoker-Phong form of the Calogero curve [24] (see also [25]) 6 



(: °*) J2 * N ~ k h% M = e '*^l f . d P Pn(X') = 0. 



.(*) 



fc=0 



,{z) 



(68) 



The Calogero-Moser-KP "quasimomentum" (that is the meromorphic differential with single first order pole 
and vanishing A periods) is 

dp = d\' = d(\- %¥£) (69) 



e*(z) 



and the Seiberg-Wittcn periods are 



a ; = <p log \dz, a i 



log Adz, 



B, 



in the Ruijsenaars model. The period matrix is, as always 



13 da, ' 



i,j = l,...,N. 



(70) 



(71) 



2.3 Theta Functions for the Calogero-Moser Family 

We will next consider several simplifications that arise when the spectral curve (12) covers a 1-dimensional 
complex torus, or its degenerations, a cylinder or sphere with two punctures. This setting contains the elliptic 
Calogero-Moser and Ruijsenaars-Schneider models, corresponding to SW theories with adjoint matter, with the 
degenerations including the periodic Toda chain corresponding to pure gluodynamics. In the case of a curve 
covering a torus the theta functions defined on the Jacobian of the curve simplify, and these are the expressions 
we shall need. 

The basic Riemann theta-function function with characteristics is defined by 



e 



(z|t) = y e 27ri s™i(™.+^-)(2.+4)+^E,™=i(™.+^) T . 3 K+¥) 



(72) 



*eZ' 



and we set 9 








(z|T) = 0(z|T). By taking T to be the period matrix of a curve of genus g = N we associate 
a theta function to the Jacobian of the curve. The period matrices of our models satisfy the constraint (22) 

N N 






(73) 



3 = 1 



and this allows simplifications. We will show that 



N-l 



9(z|T)=2^(«|iVr)e fc (4|f) 



(74) 



k=0 



6 In (64) we have chosen cr-functions for convenience instead of (?*-functions, this differs slightly from the general formula (68). 
However, one may easily check that this difference is inessential and the result is still 



Jl_ 
OX' 



e,(z) 



(A' 2 - A'/if M + h° M + 2 V ) + 0(e 3 ) = e 2 (\ 2 - Aftf M + h° M - p{zj) + 0(e 3 ). 



(67) 
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where 6±{z\Nt) = 9 



2k 

N 





(z\Nt) is a genus g = 1 or Jacobi theta-function with characteristic, 0^ is a genus 



N — 1 theta function and z, z and T will be defined shortly. 

Let ejy = (1, ... 1) be the A 
Because of (73) one may write 



Let e^v = (!,...!) be the N vector with all l's, and from this construct the projection matrix P — -h^e 



T = tP + T, f= (l-P)T(l-P). 



,jV rfi V J n V» v^AT 



Thus Tij — jj + Tij and X^j=i -^ij = = X^i=i ^j- We may similarly decompose 

z = -e JV + z, z = (l-P)z. 



(75) 
(76) 



^A" 



,7V - 



Thus Zi = jj + Si and YJ i=1 z « — z or Si=i ^« = 0- I n order to express 0(z|T) in the form (74) introduce the 
matrices 

/ 1 ... -1 \ 

1 ... -1 

M= : •.. : 

1-1 

\l 1 ... 1 1/ 

and the change of basis 



M~ 1 = - 

N 



( N-l -1 
-1 N-l 

-1 -1 



-1 1 \ 

-1 1 

N-l 1 

-1 1/ 



Mz, 



where z is now an N — 1 vector. Then 



N-l 



n.z = n J M" 1 Mz = > (n,- )zi H = (n ejv_i).z H 



j=i 



A 



with A: = y^ ■_ 1 rij. Also 



iV 



AT 



n T Tn = r^M" 1 M(tP + f)M T M T ^n = — k 2 + n T ' M' 1 MTM T M T _1 n. 



Now 



MTM T = M(l - P)T(1 - P)M T = 



T 







V ... / 



where (for i,j<N—l) 

Tij — T^ — T iN — T N j + T NN — T,^ — T iN — T N j + T NN = T t j + 2_^ {Tik + T k j ) 
Thus 



N-l 



fe=i 



Af-l 



E^ 



k,l=l 



n T Tn = ^k 2 + (n - - ejv _i) T f (fi - -e^i). 



In terms of these quantities we see 



„27rm.z+27rn Tn 



E 



0(z|T)= JT 

^+» W ^fc 2 +2«i;^ri 1 (ni-w)^+^E|^ii("l-w) f im(n m -^) 



»6Z W ii:£ini=fc 



JV iT '" Af 






T E ■ 



27ri(n-iejv-i)z+iTr(n-ieiv-i) T(n-iejv-i) 



(77) 



(78) 



(79) 



(80) 



(81) 
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By writing k = Nm + I with me 2 an d I € 2w = Z modiV (i.e. i = 0, 1, . . . , TV — 1), one has 



0(z|T) = ^ ^ e 2 7 ri(m+-i r )^+i7rA r T(m+ 1 ^) y^ e 27ri(n- -^e N _i)z+i7r(n--^ r e]v-i)' r f(n- iejv-i) 



^Z N ^Z 



neZ' 



That is 

6(z|T)= ^ ^(«|JVT)e,(a|T) 

where #j_ is the genus g = 1 theta-function introduced earlier while 



(82) 



6;(z|T) = ^ e 27ri(n--i ; -e J v-i)i+i7r(n--^e i v-i) r f(n--i-e A 



e 



ieZ ? 



2/ 



ejv-i 







(-z|T) (83) 



is defined on a (TV — l)-dimensional complex torus. (For N — 1 > 4 this torus corresponding to T, is not 
necessarily a Jacobian.) Our expression (74) is equivalent to that obtained in [23, 26] upon observing 



e fe (z|f) = Y^ e 2 ™Z?=i n i*o+™T,L'=i n J f ii' n i 

neZ™;i:f =1 ™ 3 =fc 



imj T jjimj , =e fe (z|f). 



(84) 



We have therefore established (74) for those models with period matrices satisfying (22). The coefficients 8^ 
will appear in our later discussions. 

A B 



For future reference we note that under the transformation T = 



C D 



G Sp(2g, R), where A l D - 



C B = l g , AC — C A and B D = D B, we have the arguments of the theta function transforming as 



T -» T r = (AT + B)(CT + D)- 1 
z -» z r = [(CT + D)- 1 ]*;? . 



(85) 



Further, the characteristics transform as 



e r = L>e - Ce' + ±diag(CD*) , 
e' r = -Be + Ae' + ±diag(AB*) 



(86) 



r r , M 

ForT = I m -it I \«' s< ' (> 



T r = MTM T = 



T 



° \ 





z r = / z 

z 



V ... Nt J 



and our previous discussion has simply used this transformation to put the theta function into a canonical form. 
Further the particular case 



9(z|T) = C(detT) 



-1/2 -i7TZ T T- 1 Z Q 



<d{T- l z\-T- 1 ) 



(87) 



(for some £ 8 = 1) leads to 



e(z|T) -E^o 1 © 





-1/2 



(z|JVt)9 







(-z|T) 



= C dct T (Nt) 



-1/2 i „zf- 1 z.-i-K4 r - v^iV-l o. 



EL7 @(^ + #|#i) e(r-^ - a^i - r- 



(88) 
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2.4 Elliptic solutions 

We conclude this section by examining some explicit solutions to the intcgrable systems described earlier. 
The simplest non degenerate periodic solutions arise in integrable systems with only two interacting particles. 
Because the center of mass decouples in the Calogero-Moser families we have been considering there is effectively 
one degree of freedom in this case. We will consider the "periodic" (sine-Gordon) Toda chain and Calogero- 
Moser models. 

For the example of the "periodic" Toda chain with two particles an explicit solution is a simple consequence 
of the addition formula for the (Weierstrass) elliptic functions 

p( M + + p(m - - 2p(/i) = -|^ log {p{n) - p(0) • (89) 

Indeed, consider 

e-^ + ^e-- (90) 

where fi — Ut, uj is any of the (half-) periods ( p(/j, + 2u>) = p(/u) an< 4 f] = C( w )) an d the constants A and U 
have yet to be determined. Then 



e 



'' = ^^fee^ = -# *%$X? = A2 (PM - **"» ' (91) 



^ 2q = ^" 2 pM - pM - urn* (p(m + «) - p(w)) 



Here i/ 2 = (e — e+)(e — e_), with e = p(u>) and e± = p(w±), where w± are the remaining half-periods: 
u> + lo + + uj- = modulo the lattice. If one puts A 4 H 2 = 1 then 



e 



:1 " - <■- ' 2 " = Tj (p(M + w) " Pip)) = m (P(M + u; ) + P(J* - u )- MaO) = -53? sj? lo § (*>(/*) - &>M) 

_ _J_ a 2 !»„ o-(M+")o-("-M) _ _J_ S 2 l no . p 2,j _ _J_ 9 2 <? 1 d 2 q 

~ 2H a^ 1U & <T 2 (/j)<r a (w) ~~ 2H a^ 1U «= e ~~ h a^ 7 ~~ UTT 1 IF ■ 

(92) 
After introducing A via 

u -dn ,93) 

formula (92) acquires the "canonical" form 

^| = A 2 (e 2 <? - e" 2 «) (94) 

with A 2 playing the role of a coupling constant. We may rescale A — ► 1, in which case these are the equations 
of motion for 

h=;p 2 + e 2 « + e- 2 «=p 2 + 2cosh2<7 (95) 

which is the periodic Toda Hamiltonian (33) in the centre of mass frame. (Here q = q\ = — q^.) The Lax 
operator for this example is 



C{w) = w j (96) 

and the corresponding spectral curve in this case 



e q + we 9 —p 



w+-=\ 2 -h. (97) 

w 

The scale parameter A may be restored in the above by setting p — > p/A, A — ► A/ A and w — ► w/A 2 . Finally the 
perturbative limit leading to (6) is obtained by setting p — ► p/A, A — ► A/A, w — ► w/A 2 and also q — ► <j — log A 
and taking the A — > limit. Another convenient representation for the spectral curve (97) is the elliptic 
parameterization 

w 3 + hw 2 +w = \ 2 w 2 , (98) 

where 

ct(z + w)cr(z - w) 
w = x-e = p{z)-e = p{z)-p{u>) = — 

2 , , , A (99) 



a(z + u) c _ c(u)z \ = / a(z - cj) cC(u) 
a(z)a(ui) J \<7(z)(t(uj) 
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and 

A=i^i£L. (100) 

2 p(z) -e K > 

In this form h = e, e + e_ = 1 and g-2 = 4(/i 2 — 1), 93 = Ah. Also dz = 2-^j. Observe that A — ► 00 as w — ► 0, 00 
(equivalently z -> w, or 1 ^ e, 00). Let us denote P + to be the point A — > 00, w — > 00 and P_ to be the point 
A — > 00, to — ► 0. Further let P(l) be the point A = p, to = — e~ 2q and P(0) the point A = —p, w = —e 2q . 

The solution we have just obtained may be derived from the Baker- Akhiezer function for this problem. The 
auxiliary linear problem C^ = A$ has solution 

\ / 1 \ 

(101) 



while the "conjugate" equation ^>C 



Consider the expression 



*=(::)= 


( ' 

| we q + e~ q 


\ A+p 


A\& has solution 




* = (^o,V'i) = 


fwe q +e~ q 


{ x-p ' 


ipi we q + e~ q 


V'o 


A+p 



(102) 



A 2 as A — > 00 it has a pole at P + . Also there is 
a further zero at P(l) and a pole at P(0), and so we have the divisor of ipk 

(V» fe ) = P(k) - P(0) - fcP+ + fcP_ . (103) 

Now a generic Baker- Akhiezer function should have the same number of poles and zeros. Introducing 

Vfc = U7 fc/ Vk ( 104 ) 

gives 

^=«,V>i (105) 

V'o Vo 

and 

(^i) = P(») - P(0). (106) 

In particular we may express the the Baker- Akhiezer function in terms of the elliptic parameterisation (99) as 

~ 6 (z - z(P((i)))) , s 

^ = const - o\z-z(pmy (107) 

We conclude this section with the simplest example of the elliptic Calogcro-Moser model, that corresponding 
to two particles. In the centre of mass frame /if M = p\ + P2 = eqn. (42) turns into 

A 2 + h - x = 0, (108) 

where h = —h^ M — —p\P2 + p(<Zi2), x — p(z) and we have removed the m 2 dependence in (42) using the scaling 
properties (see Appendix A) of the p function. This equation says that to any value of x one associates two 
points of S CM 

A = ±Vh=x, (109) 

i.e. it describes E CM as a double covering of an elliptic curve ramified at the points x = h and x = 00. In 
fact, x — h corresponds to a pair of points distinguished by the sign of y, but x — 00 is one of the branch 
points, and so the two cuts between x = h and x — 00 on each sheet become effectively a single cut between 
(h, +) and (h, — ). The spectral curve E C ' M may therefore be considered as two tori glued along one cut, and 
so £^f 2 has genus 2. This is a hyperelliptic curve (for N = 2 only!). The two holomorphic 1-diffcrentials on 
£ CM {9 = N = 2 ) can be chosen to be 

dx \d\ dz dx d\ 

dv + = dz = —— = GM7_ = — - = — — = — 110 

+ 2y y X 2yX y K > 
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so that 

dS = 2Xdz = A— = — \Jh - x, (111) 



\dz = 


dx 

y 


y 


- X. 


ddS 
dh 


^, dx 
= 2yA : 


= dv-. 





and 

r)d.S dx 

(112) 

The fact that only one of the two holomorphic 1-differentials (110) appears on the right hand side of (112) 

is related to their different parities with respect to the ^2 ® Z2 symmetry of H CM : y — ► —y, A — ► —A and 

du± — > ±di>±. Since diS has a positive parity, its integrals along two of the four elementary non-contractable 

cycles on S automatically vanish, leaving only two non- vanishing quantities a and a , as necessary for the 

SW interpretation. Moreover, the two remaining nonzero periods can be defined in terms of a "reduced" curve 

of genus g = 1 

3 

Y 2 = (yX) 2 =A(h-x)\{{x-e a ), (113) 

a=l 

equipped with dS = (h — x) %■■ This curve arises when directly integrating the equations of motion since after 
decoupling the free motion of the center of mass we again have a dynamical system with only one degree of 
freedom. From the conserved energy h = p 2 + p(q) wc obtain 

/da f dx 

= = / - (114) 

\Jh - p{q) J yjA(x - ei){x - e 2 ){x - e 3 )(h - x) 

which is exactly the Abel map of the reduced curve (113). 

3 The weak coupling limit 

In this section we are going to study the weak coupling limit in some detail. First we consider the A —>■ limit of 
the periodic Toda chain, demonstrating that this leads exactly to the open Toda chain or Toda molecule. Here 
we explicitly derive formulae for the theta-functions (tau-functions) in this limit and show that they appear to 
be a finite-dimensional analogue of those appearing in matrix models. We discuss their relation with duality 
in integrable systems and the commutativity of theta-functions [26] . We then turn briefly to the corresponding 
properties of the weak coupling limit of Calogero-Moser models. 

3.1 From the periodic Toda chain to the Toda molecule 

The perturbative limit corresponds to r — > +ioo and A — > 0. The effect of this will be to simplify the theta 
function solutions describing the general finite gap situation. 

First, in the limit r — ► +ioo, the Jacobi theta-functions in (74) turn into exponentials and one finds that 

e(z|T)= J2 ^ iJkz ®k- (115) 

We may already deduce an interesting consequence: the ratios of the coefficients Ofc = 0fc(z|T) = 0fc(z|T) 
Poisson commute, 

|i,|l} = „, y ,,.,,,,,. (11 „ 

The Poisson bracket here is that corresponding to the symplectic form 

iV-l N-l 

n Toda = Y^dz t A d ai =J2 d ^ A d P*> ( 117 ) 

i=\ i=l 

where q"i and Pi are co-ordinates and momenta of the Toda chain and we are working in the centre of mass 
frame q\ = qi — qcM- The Poisson commutativity of the ratios (116) follows from the solution of the periodic 
Toda chain [27, 28] 

e* = ^-, |i= J! e "- (118) 

k=j+i 
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(Here O = <9n and the @k can be thought of as the Toda chain tau-functions depending on the discrete time 
fc, the number of the particle. Observe J^iQi = here.) Now because the coordinates qi obviously Poisson 
commute, {qi, cjj} = 0, we deduce (116). This expression gives a precise formulation of the old expectation that 
the Toda chain tau-functions Poisson commute with each other. We will return to this point shortly when we 
discuss duality. Henceforth we shall drop the tilde from the coordinates and momenta for simplicity. 

The second set of simplifications arise because the spectral curve (38) degenerates into a rational curve (1) 

N 



= PN(\) = H(\-a t ) (119) 



(X) a i ~ 0) hi the perturbative limit A — > 0. Then upon making the natural choice (for the cylinder) z = log 



«' 



'At w 

and we have the basis of (normalised) holomorphic differentials 



A — (120) 



dX 

2niduj i = . (121) 

X — ai 

d"-? /\„U™„ „D _ X \ dw 



The perturbative period matrix Tjj = -^ ± - (where a® = § B A—) is given by 

-infij = -iirfrr = % J>g J^i - (1 - 6a) log M. (122) 

We shall substitute this into (115) and take the A — ► limit, but first we should be more careful with the 
various appearances of A. We have already mentioned the scalings hk — * hk/[i k , w — ► w//j, n along with 
A — » X/ ii. Comparison of (115) and the relation between w ~ e z+47rr and z shows the first two of these scalings 
to be achieved by 2itiz — > 2itiz — Nlogfi. Now the double scaling limit (where A = me' l7TT / N ) is achieved by 
shifting 2iriz — > 2iriz — N 2 \ogA. (We shall justify this a little later.) Upon noting the r dependence of the 
genus one theta function we see &k — * A fe 0& and overall we must construct 

9 fe = lim A k2 - Nk @ k = 
EiLi"i=fe; Ei<j(«*-«j) 2 =fe(JV-fe) ( 123 ) 

_ V^ e 2-ni'£™ =1 n j z j -i-!r'£™ <j = 1 (ni-n j ) 2 log\a ij \^ 

The quadratic constraint here appears when the A dependence of Tij is taken into account. Now the conditions 
"^rii — k and X^<7'( n 2 ~ n i) 2 = k(N — k) can only be satisfied for each m £ {0, 1}, with k of these nonzero. 
Thus we may write 0fc as 

e-En^no d24) 

|/|=fcie/ ,-e/ 1 ^ 1 

for some set of indices I = {i\, . . . ,ik}, and / is the set of indices complementary to / (m = 1 if i G / and 
rii = otherwise) . These are expressions for the tau-functions of the open Toda chain or the Toda molecule. 

Let us consider the N = 2 version of these formulae. For this case of one degree of freedom we have the 
hamiltonian (6) H = p 2 + e 2q . Taking for this example the period matrix 



/ , K -a 2 \ |ai -a 2 | \ / , a a 

Off — Off 



-iirWTij] 



lo>5 — A — ~ l0g A 

V -log^^ log^i^ y 



log - - log - 



A °A 

a a 

A l0g A 



(125) 



formula (123) gives in this case Oo = 1 and 

=( i -j) z - i7r ( i -i) 2 i°g a 
»+j=i;(»-j) 2 =i 



©i= E ' 



(126) 
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where we have substituted 2irizi — ► z, 2iriz2 — ► —z. Thus 



6! 



Y^ e -(i-]) 2 \°g a+(i-j)z _ }_iz _|_ e -z\ 



i+j=i;(i-j) 2 = i 



(127) 



Now e 2q = e 92 9l = (62/61 ) 2 = 1/0? and so we again arrive at the explicit solution to the equation of motion 7 
for the hamiltonian h = p 2 + e 2q 



coshz 



In the N = 3 case eq. (124) gives 



(128) 



27rizi 



01 =e 



and one further nontrivial function 



012013 



+ e 



2772^2 



|ai2a23| 



2-7riz3 



|ai3fl23 



(129) 



Q _ e 27ri(Ji+z 2 ) 1 e 2iri(zi+z 3 ) , e 2iri(5 2 +Z3) 

Iai3a23| Iai2a23| I012O13I 



_ -27ri5i _ 



1 



+ e 



— 2X222 



1 



— 2niZ3 



1 



(130) 



|ai20i3| |ai2a23| |ai3023| 

There is a convenient determinantal form for formula (124). In general one has that 

fe = det K n+m 

kxk n,m=l,...,k 

where the "moment matrix" K nm = K n+m is defined as an average with respect to ©i, i.e. 

N N 

la.j 



K n = {a n ) 1 = J2 eZi(X 



n-2 



e z > = e 2 



i=i 






(131) 



(132) 



i¥^i 



?=i 



The proof is similar to that in matrix models [29]. Indeed, 



0* 



z2\i\=ki[iei e niZl Tlj£i |o y 1 - Z_,|/|=fcILe/ e * "* [Iljei \a zj \iljei\{i} \a zj \Ti-jei\{i 



y\Uij 



— i2\i\=kiiiei ( e i rijG/\{i}i a yi ) — z_/|/i = ferii e / e *rii^ie/i a yi ~ 5_/i/i=fcrii e / e *ni<,-; ;.,•<=/ 



,2iTl n 2 



(133) 
looks like a "discrete" analogue of the tau-function of the "forced" Toda chain hierarchy which plays a central 
role in matrix models [29] . Now 8 for any coefficients Ci 



Z^\l\=ki.iiel Q lli<j ; i,jel a ij — t-il: h<...<i k ^U • ■ ■ ^*fc lli„<i m \ a i^ fl > m ) 

= E 4l <...< lfc C»i • • • C» fc det fexfc a^ 1 | , eJ . n=1 _ fe det fexfe af | . £J; m=1> _ fe 
- dct fcxfc feili C l0 Jl +m - 2 ) 

V / n,m— l,...,/c 

Substituting Cj = e z % one arrives at eqs. (131), (132) 



(135) 



I: ii<...<ik 






= det if„ +m . 

fcxfc 



(136) 



=i,...,fc 



These formulae in fact coincide with the solution of the TV-particle open Toda chain problem in terms of the 
action-angle variables discussed in [17]. 

7 It is interesting to point out that the equation of motion for the co-ordinate Q = a = \/h in the dual system H = co °r 
coincides with the equation of motion for the rational Calogero model, Q = -^j. Similar dualities are known for the Calogcro 
models and their "relativistic" counterparts. 

8 This is a particular case of the Cauchy-Binct formula for the k X N (k < N) rectangular matrices A n i and Bi m (i = 1, . . . , N, 
n,m = 1, . . . , k) 



det y^ A ni Bi,, 

kxk \ *— i 



1=1,...,*; *l<---<»k 



Edet A ni 
kxk 



det Bi„ 



itzl] n=l,. .. ,k 



i£I; ro=l, .. ..k 



(134) 
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Let us relate these formulae to the Baker- Akhiezer function. First, we have defined the angle variables as 
co-ordinates of the Jacobian of genus g = N curve E* by (19). For the holomorphic differentials (121) we have 

N rP* dX ~ " rPu dX 



Here Afc = A(Pfc) and 



, P A -a, ^ l f-^J Po A" +1 A fc ^A(p ) ., 

fc=l ° n=0 fe=l - r ° n=l 



JV 

n 



1 JV 

T « = --Ev ( 138 ) 



fc=i 

is the so called Miwa parameterization of the "canonical" Toda times. The logarithmically divergent first term 
here is absorbed into the renormalisation z% — > z% — j^logA. Summing over i yields z — > z — j-^logA and 
the double scaling limit we have already mentioned. The remaining "finite" part of the right hand side here 
describes the Toda chain dynamics with respect to the various higher times t n , n > 1 

2irizi = a i t + '^2 0?*n + zf' . / 13g \ 

n>l 

The ordinary time t\ = t here corresponds to the evolution with respect to the Hamiltonian quadratic in the 
momenta (or quadratic in canonical action variables) 

2mz, = a,t + zf ] . (140) 

The Baker-Akhiezer functions on the degenerate spectral curve (119) may now be defined by the following 
analytic requirements: It is the set of functions ipk — ipk{X), k — 0, . . . , N — 1 which have exactly fc zeroes on 
the rational curve (119) and a single pole of order fc at A = oo. For k > N they can be defined by ipk+N = wpk- 
This means each ipk(X), k = 0, . . . , N — 1 may be represented by a polynomial with k (finite) zeroes. Thus ipk 
can be constructed as a linear combination of J\ i < <i (A — Oj,) = n,e/(^ — a j)- ^ n ^ ac ^ one nas 

, m ^ fc ( Z '"^^>^) Sm =fc n, e /e M "(A-a,)n je /T^T mn 

^ kW = eTT^ = — v ff — P 2^n - * ' ( ' 

Thus, for example 



V-i(A) = A- 



6 'W siLi^iw 



i#* faTJl 



= A- 



Ef = i«,e 2 — n 



(142) 



■vAT 






V JV p27riz;TT 1 

2^=1 e lljjti]^J\ 

with (a) = (a 3 ) 1 /(a 2 ) 1 . We note that these Baker-Akhiezer functions satisfy 



A ! 



which is a consequence of the identity 



J2^Ma l ) = 



N 1 



a, — a 
1=1 j#« 



^JV 



(This may be established using Lagrange interpolation: 1 = Ei=i Y[j=u ._ s . •) One may also easily write the 
equations of motion for the zeroes of the Baker-Akhiezer functions ipk(%) — 0: 9 

^ = nU^- a i) ** = ir'uU^-^ (144) 

^ IIj^(7i-7?)' ^ Il^ifri-Tj) 



9 They can be considered, for example, as a degeneration of corresponding equations for the zeroes of the Baker-Akhiezer function 
of periodic Toda chain 



9 * n^»(7»-7j) : 9*n n^i(7i-7j) 



(143) 
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Finally let us note that the conserved Hamiltonians are straightforwardly given in terms of the minors of the 
Lax matrix. The Lax representation gives 



N N 

a,i 



det(A -C) = JJ(A -cn)=J2 X N - k (-l) k h k 



with ho = 1, and 



h k = ^2 ai x ...<Xi k k = l,...,N. 

ii<...<i k 

On the other hand, since (for any matrix C) 



det (A - C) 



X N _ X N-1 £ £ .. + x N-2 £ det ( ^ ^ \ x N-3 £ det 
i i<j V *-Ji *~33 ) i<j<k 



f^ii f^ij f^ik \ 



A' 



= \ N - X n - 1 TtC + ^(-l) fe A Ar - fe J2 det 



fc-2 



i-ji .Ljj i-jfe 
\ £fci £fcj -Cfcfc / 

rr r 

\ L-ikii >~'ikil ■ ■ ■ *~ikik / 



we have that 



h k = ^2 a h ■ ■ ■ a i k — X! dct 



' A-siiii *^iii 



C C 

*"22«1 *~22«2 



\ *-i*»i *-» 



t-ikik / 



(145) 



(146) 



(147) 



(148) 



3.2 Theta-functions for the Calogero-Moser models and their Poisson Commuta- 
tivity 

In the preceding discussion of the Toda chain we saw that the ratios 10 of the theta functions <d k Poisson 
commuted by appealing to the explicit form of the equations of motion. In fact a more general result holds 
for the Calogero-Moser family that we shall now describe. We have seen that the theta functions for the 
Calogero-Moser family satisfy (74). We wish to show that for appropriate solutions (and for all k,l,m, n) 







Ofe 9„ 

e? a; 



9/ Skr, 



9fc Sir, 



(149) 



where 



Simn = B;{9 m , 9„} + 9 m {9„, 9;} + 9„{9z, 9 m } = S mn i — —Si nm - 
Indeed (upon setting I = n here and using the antisymmetry of S) we see that (for all k, m, n) Skmn = or 



N 

o iJkmn — / y 
a=\ 



9/c 9 m 9„ 

ae k de m de n 

dq a dq a dq a 

dO k d<3 m ae„ 

dp a dp a dp a 



Such addition formulae are very restrictive [30] and closely connected with integrable systems. 

We establish for the Calogero-Moser system the commutativity (149) in the following way [26]. According 
to [19] the equation 

= 9(z|T)= J2 °^ z \Nr)Q t 



^TL t 



(150) 



10 It is perhaps helpful to recall at this point that theta-functions are used to embed a curve into some projective space giving 
the inhomogeneous co-ordinates of the embedding. Their ratios may be considered as "normal" or homogeneous co-ordinates. 
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(as an equation on the z-torus) has exactly N zeroes jf = qi, ■ ■ ■ ,qk- As a consequence one gets a system of 
linear equations 

JV 

J2^(N qj \NT)Qi = j = l,...,N. (151) 

i=i 

The system should have nontrivial solutions, i.e. dety 9_L(Nqj\Nr) = 0, which effectively reduces the number 
of degrees of freedom from N to N — 1. Then (151) can be rewritten as 

J20^(N qj \NT)^=e t o(N qj \NT), Vz ; j = l,...,N-l, (152) 

and so using Cramers rule 



0. _ det k jiio,i-^i ;j=i,...,N-idk r (Nq j \NT) 
Q io detk 7 a ;j=i,...,N-iOk(Nqj\NT) 



(153) 



Therefore, the ratios H 1 depend only on the co-ordinate qk, k — 1, . . . ,N oi the Calogero-Moser particles and 
so obviously Poisson commute with respect to the Calogero-Moser symplectic structure 

N N 

qCm = ^2 dq . A dp . ^^dz.A da, (154) 

restricted to Ylj=i Qj = cons t- The latter condition comes from the of vanishing the determinant 

det 6 i (Nqj\NT) = 0. 

ij N 

Indeed, using the ^-function identities described earlier coming from Wick's theorem [22, 21] 

det0£(Mfc|M-)~0 E £ iN^qklNr) Y[$.(N qi - N qj \N T ) (155) 

v V k J kj 

and we can can compute (153) explicitly; the vanishing of the determinant corresponds to X) 7 =i Qj being a zero 
of the theta function # s j_ (with characteristic being the sum of the characteristics jjr). This centre of mass 
constraint is also equivalent to (23), X^=i a j ~ const. A consequence of the the ratios S 1 depending only on the 
co-ordinates of the integrable system is that they may be used to construct a set of independent hamiltonians 
for a dual system [23, 26]. 

3.3 The Perturbative Limit of the Calogero-Moser Models 

The elliptic Calogero-Moser model degenerates in the perturbative limit of the SW theory r — ► +ioo giving rise 
to the well-known trigonometric Calogcro-Moscr-Suthcrland model. The solution in terms of the action-angle 
variables is a direct generalisation of the open Toda chain case and may be presented in terms of the dual 
rational Ruijsenaars-Schneider model. The salient features arc as follows. 

The Lax operator 

£ij = PiSij + m(l - 5ij)-— — (156) 

Slllil Qij 

may be considered as a limiting case of the Lax operator of the trigonometric Ruijsenaars-Schneider model (50) 
in the same way as we derived the Lax operator for the elliptic Calogero-Moser model from its Ruijsenaars- 
Schneider counterpart. The spectral curve for the model is a minor modification of (119). Indeed, from (68), in 
the limit r — > oo one gets (for appropriate imaginary period) 

^f-"^W ) = 0. (157) 

smn(iz) 

Upon introducing w — e z , one may express this as (8) 



w 



Pn{\) A A- a 



Pn (A — m) ■*■ *• X — a 

2 — 1 
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-*- J- A — n ■ — m 



Again for su(N) we have J^ CLi = 0. The spectral curve is equipped with a generating differential dS = A—. 
The period matrix in this case is given by [24] 



,2 , a 2 



_ inT Per, = _ iwrgjk + l tf , fc £ lQg (Qrj+ro °^_ m) - 1(1 - **) log (o . fc+Tn ^ Q . fc _ ro) (159) 

where the first term corresponds to the bare coupling of the elliptic Calogcro-Moser model, the modulus of the 
base torus r = 9 -\ — ^ — . In the perturbative limit r — ► +ioo (159) this is renormalised and remains finite. In 



9 



YM 



N 

2irikz 



this limit one gets from (82) 

6(z|T) = J2 e^e^&k (160) 

fc=0 

so that, finally, the "tau- functions" of the trigonometric Calogero models may be introduced as 



|j|=fc*ei -jei 
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4 The Strong coupling limit 

4.1 Solitonic solutions of the periodic Toda chain 

The quantum moduli space of the 4D pure su(N) J\f = 2 SYM has N maximally singular points at which TV — 1 
monopoles become simultaneously massless. These are the confining vacua of an Af = 1 theory [1, 31, 32]. 
At these points the dual variables af are the appropriate variables to describe the prepotential. The Af = 2 
spectral curve (32) at these points is described in terms of a Chebyshev polynomial P N c ys CV (A) defined by 
Pn ° yS ev (2cost>) = 2cos(Nv). With w = e z and this choice of polynomial we see (32) turns into 

P^ hcbyshcv (A) = 2coshz. (162) 

From the definition of the Chebyshev polynomial it is clear A = 2 cosh (-^) is a solution of (162) as indeed is 
A = 2 cosh (-^ + i^f-) (k = 0, . . . N — 1). These are the Af = 1 points of the theory, related by a 2w symmetry, 
and we will focus on the first of these in performing our analysis. 

The hyperelliptic form (35) of the spectral curve (recall y — w — l/w) is now 

y 2 = P N (A) 2 -4=(A 2 -4)Q(A) 2 (163) 

where the roots of polynomial Q(X) are given by 

JV-l 

Q(A)= ]7(A-2cos^). (l 64 ) 

3=1 

This is a "solitonic" curve in the periodic Toda chain: if we express the curve as y 2 — Y\j=i (^ — e fe) we see that 
e2k = £2k+i = cosirk/N (k = 1, . . . N — 1) and the corresponding B periods have collapsed; e\ = 2 = — e2N+2 
are single branch points. 

Let us now introduce a new variable by 

A = 2cosh-^ =£ + <T 1 . (165) 



Now (165) maps the 2-sheeted cover of the A-planc y = VA 2 — 4 to a cylinder with co-ordinate £. Thus eqs. (162), 
(163) describe analytically a cylinder with N — 1 distinguished pairs of points. With these coordinates our 
differentials on the curve now take the form 



x^d\_ x^d\ m N - 2 ^ + l) k ~ l 



y VA 2 



4 nf = - 1 1 (A-2cosf) qJL-^-eV) ^ 



— e N 



,7V- 1. (166) 
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For a non-degenerate curve these were holomorphic but now they acquire simple poles at the points £J~ = e « 

and £~ = e~~«" (and have singularities at the infinities of the cylinder (165) £ = 0, oo). The canonical basis in 
the space of differentials (166) can be chosen as (j = 1, . . . , N — 1) 



D _ sin ¥ d£ _ sin w d £ 



^ " ( e _ e ^)( e _ e -^) tt (£-#)(£-§-) 

These differentials are normalised to the B-cycles, here the cycles around the marked points O 

,d _ I j, n 



(167) 



«+ HI 



dwf = <b duj = - f duij = 6ij, (168) 



while certain of the A-periods (§ A du® — J i ckvf) diverge, the others (j ^ fc) being given by 

3 So 



T ,k = f duj k =7T- log - 2 Z L\. , ,, = — log . t , ■ | ,. v (169) 



J_ sin 2 ^(j -fc) = 1 sin^lj-fcl 

2™ ° S sin 2 ^(j + k) m ° S sin ^ (j + fc) ' 

Using this expression one may show that T D (169) satisfies the Edelstein-Mas [10] identity 

Erfi' Trkj' v-^ ~nert, 7r(' — I) n 7rfc(« — 7j) irkd — h) N 2 ^ n 

sin sin — rf- > 71?- (oj -» 2 cos — ^) cos — ^- cos — ^- = — TR, . (170) 

k—1 «J— 1 

This conjecture of Edelstein and Mas is proven in Appendix B. We note that an equivalent expression to (169) 
was also estabilished in [33] where an interesting investigation of the Af = 1 degenerations of the "multisoliton" 
solutions is presented. 

The Abel map in the present setting is 



JV-i ,?, , iV-l 



z 



? = y rdu? = -Lviog^^-^iy;t n 8in^ 

fc=l fc=l SK S J n=l 



E ^-^E^^-iE*.*^ (»") 



where £ = {£fc}, fc = 1, . . . , TV — 1 is the set of AT — 1 points which define the Abel map 11 . The asymptotics at 
"infinity" are given by 



and 



z? 



n=\ k K n=l k K 



Tr^^--h^r^--h^H i -i)H i 



(172) 



ik 



= ^E(«/>"-«7>")E^E^E^ 

n—1 k n—1 k 



(173) 



provided 



1 £" 

f «,= -E^T> '-rXf- (174) 



^oo fe <r $-0*^ " 



11 The Abel map here depends on the genus Af — 1 of the smooth curve (32) that arises from the double scaling limit of the genus 
N elliptic Calogero-Moser curve. It is interesting to point out that in the perturbative limit the basis of differentials (121) of the 
pcrturbative curve remembers even more - that the original curve came from the elliptic Calogero-Moser model and had genus N. 
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We observe that the relation (171) with i„>2 = 0, 

zf=t lS in^ (175) 

coincides with the vacuum value of the string tension, or monopole condensate, proportional to the SUSY 
breaking parameter t\. 

Naively, in the strong coupling limit T^ — ► (and T® — > oo) the 0-functions (83) turn into 

@k f=l E ^ fZi (176) 

which is not a well-defined object. The resolution to this may be seen by considering the N — 2 case. Here one 
appears to have 

©o =E„ 1+ „ 2 =o e2 " (ni " 1+ " 2 " 2) =E„e 2Tm ( il ^ = _ z £ n e 2 ™ z = e 3 (Z\f^0) 

Z!—Z 2 = § (177) 

1 ' ■* i i ^ ni _ n2=rl odd 

and, naively 

eg= Q 1 = , 2( Z|f^0 )=2cos7rZ 

©o e 3 (z\T^o) 

This appears to have nothing in common with the desired solitonic solution given by integrating the equations 
of motion 

*= / ; dQ = i \^S_ ( 1 79) 

J Vh - e 2 " - e- 2 " h=2 J et-e-i K ' 

in the solitonic limit h — ► 2, which gives rise to 

e q = itan(t-t ) (180) 

(for an appropriate constant of integration to). However the correct answer does appear after both a shift 
Z = Z — 2 and a modular transform T — > — A = T . Then one has 

9 i ^ l( Z|f^0) = ^i22 = itan7r ^ + 0(e ^) (181) 

where Z D = I = -T D Z = 5 + X. 

T T IT 

The combination of both a shift and a modular transformation appears to work in general. Using the 
modular properties of thcta-functions we have 



9fe(z|f ) = ^ ™»-l e 2^(m-^e)z+i 7r (m-^e)f (m-&e) 



e 



-27riiei+i7T-|iefey^ „ x g27rim(z- if e)+.™f m _ ^iri^ei+iTT^ef eg, / g _ fe_^ e 



eZ ~ " c y* - Jj 



T 



= e -27ri&ez+«r-^efe-»7r(z-£fe)A(z-&fe) ( dct yA 2 @ (f-lg _ A 



A 



-T- 



Thus 



e k _ e^-^z-Aei-f- 1 ) _ e „ e^^-Hz + iej-Aei-f- 1 



fe -i 0(f -iz - *=le| - f-i) e[e](f-i (z+ i e) _^l e |_f-i 



e» 



where now 



e[g](* p -JH-f-i 

0[o](^-^e|-f-i 



z^=T- 1 (z + ie) 



(182) 



(183) 



and e = (1, . . . , 1). The effect of the shift is to keep only the leading terms in the quadratic (m+^e)T 1 (m+ie) 
in the limit T — > yielding 

9fe - *» '*<*1 (184) 



©fe-i f^o ' ^fe-i(z D ) 
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where 



e (185) 



and Z 2 = {+1,-1}. 

Using (173), (184) and identifying the off diagonal parts of T' 1 with T D one can propose formulae for the 
Baker- Akhiczcr functions 

**(*, t) ~ f *eS. «.«•-«-) ^ ( f^ )} (186) 

where 

tffc(z D (f n )) = ^ e nZ - j ^*" J " * 42 ^ J J J Ss '°^ |j - y| (187) 

nij— ± 

and 

sin * (j+j') /< r _,-\ ™i/ 2 

tf fe (z D (i„:£)) = V e "^ £ > mi+ ' E '." m ''" si1 ' ^^^^ m '^' log -4b-''i TT [2 £L ] 

m 3 =± j \^ _ ^' / 

(188) 

^) n(if^) • ■ 

Low order cases of this are 

• SL(2). N = 2, N - 1 = 1, TOj = mi = m = ±. Then 

tf fc (z°(t n ;0)= E e -^+^^|±]J = ^^e<<*-*>fc + i) + e-'<«-*>(f-«)) (189) 



m— ± 



• 



SL(3). AT = 3, N - 1 = 2, j = 1, 2, (mi, ma) = {(++), (+-), (-+), (--)}. 

£ - e 3 / V £ - e 3 / (190) 

E. 7m ^—v . 27m 

t„sm— , z 2 = > t n sin — — . 

These formulae give rise to the general form for the solitonic Baker- Akhiczcr function of the periodic Toda chain 

**«,«) -^eE^-*-')*^. (191) 

Here R(£) is a normalisation factor, independent of times, and chosen to be a polynomial of £ of degree TV — 1 
in order for (191) to have desired analytic properties, while 

iV-l N-l 

Rk(S,t) = Mt) II (£-"«(*'*)) = E r '(fc,*)^. (192) 

s=l 2=0 

These Baker- Akhiezer functions are defined for k : . . . N — 1 and extended to all fc by 

* fe+ Ar = u;*fc, w = £ JV (193) 

(equivalently Rk+N — Rk)- Now the Toda chain Lax equation 



A*„ = C„+i* n +i+p„*„ + C n *„_i, C„ = e^ (9 " «— l) , A = £+- 

implics that 

r (n)-C rl r (n- 1) = 

ri(n) - C n ri(n- 1) - p n ra(n) = 
25 



(194) 



(195) 



and so 

r (n) = C n r (n - 1) = . . . = e^«"-9°V (0) ~ e^ 9 " . (196) 

For the solitons coming from degeneration of TV-periodic Toda chain one should also impose the "gluing condi- 
tions" 

*n&)=*n(^) j = l,...,N-l. (197) 

This means that the Baker- Akhiezer function remembers that it originally came from a genus N — 1 Ricmann 
surface, and each pair of points £j, j- corresponds to a degenerate handle. Now the conditions (197) together 
with (193) entail 

£, 2 3 N = 1. (198) 

Thus we may take 

& = e^ (199) 

where the label j can be restricted to j = 1, . . . , N — 1 since 

1 7T7 

& =& + r = 2cos -KF = ^ 2JV -^' • ( 20 °) 



& N 



Eq. (197) explicitly reads 



¥M = n f " M f ;? = e *^s. »-#+*<*>, ,=i,..,^-i. (2od 

Here Z,-(i2) = log ^|^, thus if on has chosen i?(£) = IlfLi^C - 7s) in (191) then Zj(R) = ^f^ 1 log f^. 
Now (201) is a system of linear equations for the coefficients rk(n, t) of the polynomial R n (£, t): 

E sin ( ^ + ^ + 2 E f ' sin ^ - \ z A R n rfc ( n '*) = °- ( 202 ) 

fc=0 V ^ i / 

Such can be readily solved. Conditions (200) can be interpreted as values of the scalar fields at the critical 
points of the superpotential, while the soliton trajectories connect the critical points. 

4.2 On the "solitonic" limit of the elliptic Calogero-Moser system 

The solitonic limit of the elliptic Calogero-Moser system involves many open question, though it seems that 
we can proceed in an analogous manner to the Toda chain case. We will discuss here only some explicit 
computations for low N and will postpone the computation of the soliton phases or string tensions. 

For the N = 2 case the solitonic limit corresponds to h = e = p(oo) (w may be any (!) half-period) so that 
(65) turns into 

A 2 - x + h = A 2 - p(z) + p(lu) = (A - $(z, w)) (A + $(z, w)) = 0. (203) 

Thus 

c(z — lu) /-/, ,\„ a(z + lu) M,,\, 

A = ±$ (z, u>) = ± ; . . ; e c(w)z = T ) , , { e' c{uj)z 204 

<j{z)a(ui) <j(z)a(u>) 

and these functions are related to the entries of the Calogero-Moser Lax operator (40) at half-periods via 



<j(z)<t(uj) 

F(-.\z) = -4^1 = f£±%*<«* = * ( ,, w)e -««). 

<j(z)<j{u) <j(z)a(uj) 



(205) 



In this limit one also has 

y 2 = (x - e)(x - e+)(x - e-) = A 2 (A 2 -£+)(A 2 -g_), e± = e± - e. (206) 

The Seiberg-Witten differential and periods are 

dS = X— = — dx =, <fdS = 0, (fdS = 2Tn. (207) 

y ^/(x-e+)(a;-e_) Jb J a 
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Thus instead of Seduced (113) one may introduce the rational reduced curve 

Y 2 = (x-e+)(x-e_). (208) 

Indeed, a direct integration of the equation of motion (114) gives now 

f dx 

2lt = / 7 ^ 77 v< 71' ( 20Q ) 

J (x — e)^/(x — e+)(x — e_j 

a contour integral on the rational curve (208). 

Let us note that formula (203) can be obtained by considering the extrema of the Calogero-Moser hamiltonian 
(the superpotential in the SW approach). Then 

dh dh 

dh = ^-dp + ^-dq = (210) 

Op dq 

yields 

p = 0, p'(q) = 0. (211) 

The latter is satisfied by any half-period q — u>, so that 

h\ dh=0 = p(w) - e. (212) 

A similar argument holds for the general iV-particle case. Now the Lax operator (39) computed at the "sta- 
tionary" points yields p\ = . . . = p^ = and the locus equation 



C N = {qj | q 3 + qi, ^ p'(qij) = 0, i : 1 ... TV}. 



(213) 

The (closure of) this locus has a rich geometry and many questions regarding it are still unanswered. The early 
work of [34] is still one of the most detailed investigations (see also [25] ) . The (closure of) the locus has in general 
several disconnected pieces some of which are trivial copies of the torus. The latter are easily understood: any 
odd, periodic function (with period L) satisfies the identity 

/(£)+/(!)+/(!)+-+/(^)+/(^)=° 

for any odd TV and for TV even provided / (-j) = 0. Choosing p'(x) as the function f(x) with L = 2u (and 
p'(u>) = 0) one gets a solution to (213) with 

Tfe=go + -^-fc, qjk = -j-jr(j ~k). (215) 

The loeii here are simple copies of the torus; by varying the periods in this construction one gets further simple 
solutions. There are however other solutions less well understood. 

Lets consider the TV = 3 case. For TV = 3 eq. (213) is equivalent to 

p'{qn) = p'(m) = -p'(«i3). (216) 

This has solutions (together with q\2 + #23 = 913) 

gi2 = ui, 1723 = w', gi 3 = u> + J (217) 

and 

2lu 2lu Alo 2lo , . 

912 = -y, 923 = -g", gi3 = - 3 -=-" 3 - ( 218 J 

for any half-period u> = (u>,u>',w + ui') 12 . Upon substituting (218) and pi — into the Lax equation (39) one 
gets 

\ 3 + 3F+F-X- Ff+Fl =0 (220) 

with the three solutions 

A = F+-F_, A± =e ± ^F+ + e ±1 ?F-. (221) 

Here 

F± = ^^T eT§CMz - (222) 

cr(z)cr( — ) 



12 One also has the relation between (^-functions r\ = £(uj) in the half-periods 

C(w + w') = C(w) + C(w') (219) 
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5 Discussion 

In this paper we have considered the singular limits of various SW intcgrable systems that are relevant for 
the weak and strong coupling limits of the corresponding field theories. We conclude by making several com- 
ments about the relation of these relatively straightforward calculations in integrable models with corresponding 
properties of the SUSY gauge theories. 

Let us begin with the outstanding problem of why there is a correspondence at all between integrable systems 
and SW theory Part of the problem is that, from a purely four-dimensional perspective, Seiberg-Witten theory 
only sees the commuting Hamiltonians of a mechanical system. Only these quantities appear as coefficients of 
the Seiberg-Witten curve. Half of phase space is not apparent at all and the choice of mechanical system appears 
arbitrary. We believe viewing these systems from a three dimensional perspective sheds light on the matter. In 
fact one should also pay attention to the four dimensional AT = 2 gauge theories compactified on R 3 x S 1 [35]. 
Recall that the phase spaces of integrable systems play the role of moduli spaces, or spaces of vacua parameters, 
of (compactified) SUSY gauge theories in the following way. The minima of the scalar potential in the gauge 
theories with extended supcrsymmctry 



V($) = Tr^[$ 7 ,$ J ] 



(223) 
KJ 

correspond to simultaneously diagonalisable matrices [$ 7 , $ J ] = whose eigenvalues {<p{} can be thought of as 
the (complexified) momenta of some "particles" . Now for theories with a compact dimension one should also add 
"co-ordinates" {qf} corresponding to the eigenvalues of the Wilson loops § A^dx^. For N = 2 vector multiplets 
in four dimensions 13 one has a single complex scalar and compactification of one space-time dimension gives 
rise to an extra complex scalar q — q M + ij, where q m comes the set of eigenvalues of the Wilson loop § A 3 and 
7 corresponds to the (set of) 3D dual photons d a Ap = e a /3 p d p -f, a, (3 = 0, ... ,2. Thus by viewing the theory 
on R 3 x S 1 one may naturally include coordinates. The moduli space of vacua for the theories on R 3 x S 1 is a 
hyper-Kahler manifold, and, we recall, such are the phase spaces of algebraically completely integrable systems. 
Further, it was argued [35] that there is a distinguished complex structure on this moduli space independent 
of the radius of the compact S 1 . This yields the complex structure of the mechanical system. In general we 
expect the symplectic structure to arise from a careful treatment of the central charges of SUSY algebra: 

ft ^ Z BPS ~ f Tr(F jk + iF jk )$ ~ Tr f 5 A A 6$ ~ ^ % A 6 Pi . (224) 

For a M = 4 theory, we have three different choices for the symplectic structure (224), corresponding to the 
three different scalar fields in the N = 4 theory, and these are related by a "hyper-Kahler" rotation. 

Now consider the decompactification of the "3D" gauge theory on R 3 x S 1 with symplectic form (224) in 
the limit of the S 1 radius R — > oo. The dimension of the moduli space of the R 3 x S 1 theory is twice that of 
the limiting 4D theory: the "coordinates" are no longer present. The resulting 4D gauge theory is associated 
to an integrable system: the 4D moduli are, from the 3D point of view, the Poisson-commuting hamiltonians 
constructed from the 3D momenta and co-ordinates, with respect to symplectic structure (224). From this 
perspective, quantum effects in the compactified J\f =2 gauge theory turn the bare symplectic form (224) into 

y^ 5qi A Spi i-» V^ Szi A ddi, (225) 

where the SW integrals [1] 

dS (226) 

are the correct quantum variables. One may consider a, = ai($, A) as a transformation from the bare quantities 
{4>i\ to their exact quantum values {a^} playing the role of the quantum BPS masses of the effective theory. 
In the same way one should consider the transformation q, — > Zj as transformation from bare values of the 
monodromy to the exact quantum values of the effective theory 

For theories with four-dimensional J\f = 4 SUSY the effective couplings and BPS masses (i.e. the eigenvalues 
of the scalar fields {4>i} and Wilson loops) are not renormalised, since the symmetries of the theory include 

13 An Af = 2 4D vector supermultiplet in the adjoint representation consists of an M = 1 4D vector multiplct (A M ,i/>) together 
with an A/" = 1 4D scalar multiplet (</>, x). Here i/> and x are two complex Weyl spinors. If say \ acquires a nontrivial phase (231) 
under a shift along the loop in the compact direction, it becomes massive with the mass -|-. The 4D A/" = 1 vector multiplet remains 
massless, and can be represented by a 3D M = 2 supermultiplet (Actp, ^), where a = 0,1,2, q = RA3 and i/> is 3D complex 
spinor. 
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the conformal group and so there are no dimensionful scale or mass parameters. The corresponding integrable 
model is a system of free particles: the Hamiltonians are Uk = ■^Tr<f)' £ = T^2iPii where {pi} play the role of 
momenta and the co-ordinates coming from the "compact" moduli depend linearly on the angles. Breaking 
four dimensional supersymmetry down to Af — 2 (for example, by adding extra mass terms m 2 Tr$> 2 for two of 
the three scalar fields) reduces the dimension of the "scalar" moduli space down to 27V, that of one complex 
"diagonal matrix" field. Moreover, in contrast to the J\f = 4 theory, the matrices of the scalar fields and the 
monodromies become dependent upon each other, or satisfy a nontrivial commutation relation (coming from 
the vanishing of D-tcrms) 

[A, $] - mJ (227) 

in the general case (of nontrivial boundary conditions) . Here J is some matrix of "gauge-covariant" form and 
the right hand side here is linear in the parameter of the "massive deformation" [5, 6], (for m — > one comes 
back to TV = 4 theory). 

Let us consider the compactification of an J\f = 2 SUSY Yang-Mills theory with just a vector supcrmultiplct 
to M 3 x S 1 (with S 1 having radius R) in more detail. 1 If one takes all the fields to have periodic boundary 
conditions in the compact direction this would yield an TV = 4 (in the the three-dimensional sense) SUSY 
theory. If, however, one puts [36, 13] 

<t>(x + R) = e ie (j){x) (231) 

on half of the fields, the resulting theory would have only TV = 2 t/iree-dimensional SUSY (TV = 1 in the four- 
dimensional sense), i.e. the supersymmetry will be (partially) broken by the non-periodic boundary conditions. 
Now in contrast to TV = 4 SUSY in 3D, an TV — 2 supcrsymmetric theory can generate a superpotential [37]. 
In terms of the complexified variables q = q K + ij, the superpotential acquires the form [35, 38] 



W ~ eTr$ 2 + -p-r > e gi + 1 - qi + e qi - qN . (232) 




Here the first term is the "4D contribution" while the second term (the first term in the brackets) has a 3D origin; 
the final term is induced by 3+ ID instanton contributions. The simple roots of the second term are the usual 
3D instantons (4D BPS "monopoles") and give the potential of the open Toda chain: ^{ jj/^n-i ) — ^ N ~ x - 
The final term (minus the highest root) appears only in 3+1 dimensions and can be treated as a 4D instanton 
(or caloron) contribution. In the pcrturbative limit, considered in detail in this paper, one may directly see the 
nonrenormalisability of the superpotential, which means, in particular, that 

W iD = h = a 2 = p 2 + e 2q = W 3D . (233) 

Now let us discuss how the formulas obtained in this paper may be be reinterpreted from this three- 
dimensional perspective. First, we have discussed the singular limits of the SW integrable systems, when 
the spectral curves degenerate and become rational with marked points (i.e. their "smooth" genus is zero - 
eqs. (1), (4), (8), (9), (158), (163) and (208)). Degeneration means that the discriminant of the corresponding 
smooth curve vanishes as one approaches particular (boundary) points of the moduli space. In the paper above 
we have divided these singular limits into two groups: weak-coupling from the point of view of the SUSY gauge 

14 In practice this means for the theory at mass scale A = Aqcd that A >> -= corresponds to a 4D theory while A << -^ to a 
3D theory. For the couplings one has \ = -^|, i.e. the 3D theory with fixed coupling 33 corresponds to R — > 0, or #4 — > and 

so to the weak coupling limit of the 4D theory. In this limit the 4D instantons decouple since the 4D coupling 514 is small (for 

1 

2" 

_R — > and 33 = fixed) and so that the exponential terms, powers of e a i , may be neglected. In the 3D theory the integral for the 
polarization operator 

(k^kv-S^k 2 ) — +... (228) 

m 

yields a convergent integral by dimensional arguments. In the "intermediate" case of a compactified (3+l)-dimensional theory on 
M 3 X S 1 with finite S 1 radius R one can present this result in terms of a 3D theory together with a sum of Kaluza-Klcin contributions 



9l ¥ 

For R — » only the term with n = survives in the sum (229) leading back to the 3D expression. For the opposite limit R ■ 
one can define the dimcnsionlcss 4D coupling and replace the (divergent) sum by an integral 



(229) 



1 1 _ 1 y, 1 

^ ° \n. < r>..^.. „ 
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= A da A 

— =log— . (230) 



theory which yields "open" (or "trigonometric") integrable systems; and strong coupling yielding "solitonic" 
behaviour. If one deals with the gauge theory with extended supersymmetry each of these degenerate curves 
as well as the smooth ones play equivalent roles as physical vacua. By breaking the extended supersymmetry 
down to J\f = 1 (in the four-dimensional sense) however, one generates integrable dynamics distinguishing the 
second, "solitonic" , group of degenerate curves from those of the "perturbative" ones, so that only the soli- 
tonic points play the role of vacua in M = 1 theories. It is clear from the explicit form of solitonic solutions 
obtained here that the angles, corresponding to Wilson loops, are proportional to the parameters of the partial 
SUSY violation (cf. eqs. (171) and (190) with refs. [31, 32]). For both kinds of singular limits we have shown 
how explicit solutions may be straightforwardly calculated. The co-ordinates and momenta of the integrable 
systems of particles (the classical 3D moduli) are expressed through the theta-functions of the SW curves. In 
the singular limits these theta-functions degenerate into the finite sums - see eqs. (123), (124), (133), (136), 
(161), (187) - suggestive of a kind of partition function for a "discretized" matrix models. The "dual" point of 
view is to consider the theta-functions O^ as generating functions for Hamiltonians of a dual integrable systems 
with co-ordinate's ai and momenta z%. The Poisson commutativity of their ratios, implicitly proven in [26], can 
be checked by straightforward calculation in the degenerate limits. The fact that degenerate theta-functions 
appear in the form of "discretized" matrix models (133), where instead of an integral over the real line one 
has a (multiple) discrete sum over the spectrum of the Toda molecule, may be considered as a possible (though 
very speculative) sign of a dual description of SW theory in the language of quantum gravity, possibly in a 
"holographic" sense. 
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A Elliptic Functions and the Inozemtsev Limit 

In this appendix we shall review some basic definitions and properties of elliptic functions and then consider 
the Inozemtsev limit. 

The Weierstrass elliptic function p(z|w,a/) is the doubly periodic function with periods 2oj, 2a/ given by 

p(z\u;, cj>) = - + ]T ( {z + 2ujn + 2uJ , m) 2 {2ujn + 2uj , m) 2 ) ■ ( 234 ) 

n : m£ t W 

It satisfies the differential equation 

p{zf = 4(p(z) - ei )(p(z) - e 2 )(p(z) - e 3 ), (235) 



and scaling relation 
Using 

we may write 



p(tz\tu>,tu>') =t~ 2 p(z\uj,uj'). (236) 



sin nx = ttx 



11(1-5). £77^1 = ^ 



(x + n) 2 sin" ttx 



^2 , Jl 



2 l <—i_ sin TTlz + mT) 3 *—L sin ttitit v ' 

meg v ' meg 



where z = #■ . A slight rewriting of this utilising the scaling formula yields 



p(v\ilTT, Z7r) = > x 1 > 



12 2 ^ sinh 2 (i7rmT) 2 , ^ cosh(w + 2inkT) - 1 ' 



(238) 
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This is uniformly convergent for 7m(r) > 0. 

The roots of e^ of equation (235) may be expressed in terms of theta functions by 

e i =T&(^(0) 4 + ^4(0) 4 ) 

= &* nr=i(i - « 2 ") 4 (n~=i(i + 9 2 - 1 ) 8 + n~ 1(1 - s 2 - 1 ) 8 ) 

e 2 = T f-,(e 2 (O) 4 -0 4 (O) 4 ) 

= ii^ n~ 1(1 - <z 2 -) 4 (ie«n~ 1(1 + i 2n r rr =1 (i - ^r 



e 3 =-T^(^(O) 4 + 03(O) 



= -ii^ n~ 1(1 - ? 2n ) 4 (i6«n~ 1(1 + ? 2n ) 8 + n~ 1(1 + a 2 "- 1 ) 8 



where q — e t7rT and r = — . As g — ► we have 



7T 2 7T 2 w „ „, , t 7T 2 7T 2 



ei =6^ + W)' e2 ^-T2^ + 12^ 16? + °^ * = -l2^-l2^ 16 « + W)- 



Thus in this limit 



6w 2 



e2 = e 3 = e± = 



12w 2 
tt^ I tt' z 1 (240) 

12^ 2 ' 4co 2 sin 2 7rz ~ 4w 2 ^ 3 J 

1 ,. 7T 3 COS7TZ _ 7T 3 

2 P(z)= 8^dn^i = 8^ ; ' 
y 2 = .T 2 (.T-l) 



x = p(z) = - T7 ^3 + - = — (x - -) 



Consider now m (p(i>) — p(z))- The Inozemtsev limit is a double scaling limit in that the coupling constant 
to here is scaled as well as the period of the p function. Let 

m = Me-" T / 2 , v = (f)-nrT, w = e z+ivT , (241) 

then 



171 \P\ V ) P\ z )) — 2 Z^fc=-co \cosh{4,+iit(2k-l)T)-l cosh(z+2i7rfcT)-l ) 

, 77* 2M 2 coskd>~M 2 (w + ±) . 

B Proof of the Edelstein-Mas Conjecture 

We have 

-2i7rXg crt = Sij ^ l°gK - a i) 2 - (1 - &j) logK - a,) 2 . 

Let us denote by 



(242) 



(243) 



5 = zL sm — sm — ,L T £ (ot->2co8 ^ )cos ^^cos ^- (244) 

fc=0 i,j = l 

The Edelstein-Mas conjecture is (for %' ^ f) 

S W ^T°, (245) 

and we shall now establish this. We have shown that a direct calculation yields (for j ^ k) 

n ( n 1 sin 2 ^Q-fc) 1 sin^-lj-fcl 

Tf k = i du? = — log . 2 ™Y -4 = — log . ™\ J . -i. 246 

3 J A, 2« 6 s in 2 ^j + fc m 6 sm^j + fc 
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To proceed we first perform the sum over k, 

Eirki 1 nkj' irk(i — i) 7rfc(j — A) f .r-^ ^-^ irka 1 -r-^ s—^ .irka. 
sin sin cos — cos — = — > > cos( ) > > cos( , 

fe=0 a£A + fe=0 a£A_ fe=0 

where 

A+ = {*' - f ± (i - j),i' - f ±(i + j - 1 )}, A_ = {if + f ±(i- j), %' - f ±(i + j - 1 )}. 

Now 



Trfca^ _ J i[l-cos7ra] a ^ 0,27V, 



N-l 

^ ° OS( ' N ' \N a = 0,2N. 

If every a € A± is distinct from or 2N, there is cancelling between the terms in the sums. We find (taking 
account of parity) that 



^ . nki' . nkj' irk(i - \) nk(j - \) N , ^ ^ 

2_^ Sm -JT- Sm ~iy- C0S Jr COS — = y \ / J °otfl + °a,2N ~ }_^ ° a >° + ° a 

k=0 " ' \a£A + aeA 



Without loss of generality we may assume i' > j' . Then 



N N-(i'-j') N 

X X) ( Sa >° + s °<>™) T ij >It = X ^H+[i'-j') + X T u-{i'-3') 

ij — 1 aGA+ i— 1 i—i' — j' + l 

N i'-j' 

I V "* rr^pcrt , \ ^ rripert 

+ 2-*l 1 i2N-i-{i'-j') + l "+" 2-^1 1 ii'-j' + l-i 

i=N-{i'-j')+l i=l 



(247) 



In the limit a; — > 2 cos - .y we note that 



f P." 4 = f$£ N = 2*°* = f j2Ar+1 _„ (248) 

and for i =/= j 



report 1 , . / 7T(i- |) 7r(j-j) S " 

= J_ l 0g ( 16 sin 2 *<* ~ j) sin 2 ^ + f~ 1) 
2i7r 5 V 2iV 2iV 



(249) 



To proceed we now distinguish two cases depending upon whether i' — j' is even or odd. Let us take the 
even case first. Using (248) we may combine the first and third terms of the right hand side, of (247) and the 
second and last terms to give 

N N N 

X X(^° + ^)^ crt = Z^^)+X^ CT W)- ( 25 °) 

ij — l q£A+ i— 1 i—1 

The sum over A_ similarly simplifies, though we note that here there are now distinct cases to be considered 
(depending on whether i' + j' < N or not). We obtain 

N N N 

X £ (S a ,» + 5 a a N )f?r - X f ul\i> + r) + X f u-a> + i>) (251) 

i~j — 1 aGA_ i—1 i—1 

Combining our results shows (for %' ^ j') 

N N 

n X ^ I rppert , rppert report report \ 

"" g 2-~i V ii+(i'-j') »»-(»'— j') »»+(»' +j') H-(i'+j')J ~~ 

i=l 

, T N / . 4 Tr(i'-j') . 2 7r(2»+i'— /-l) • 2 ir(2»-t'+/-l) \ 

logTTl — ^^ Sm 2W Ism 2JV \ (252) 

- 16Z7T ^ l\ \ sin 4 £M sin 2 ^iW-l) sin 2 ^(2i-«'-J'-D y 

4 *'J' 
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thus proving the conjecture for this case. Observe that no terms in this product vanish with our assumption of 
i' — j' being even. 

The remaining case follows in an analogous fashion, the difference now being that Ta terms can appear when 
i' — j' is odd. Set 8 = %' — j 1 and 8 = i' + j'. Let 8 < N (a similar argument holding for 8 > N). With 8 odd 
(250) takes the form 

N Ni N/ 

£ £ (<J a , + S a , 2N )f f^ = £ ff-\ + £ fr^s +ffZ t s ±1 + f P N er U N ^ , (253) 

i,j=l a£A+ i— 1 i—1 

while (251) becomes 

N N N 

£ £ (^,o + ^)fr = £f^ (254) 

i,.7 = l a£A- i— 1 i— 1 

Combining these two expressions yields 

o _ v*2 ~/ T'-P I / T^pcrt , rr^pert /-report rr^pert 

15 ^ ^ J i'7"' + o J -S+l S+l ' J A7 i-i jw «-i ~~ J 5+1 5 + 1 — 1 



i'j' -r 8 y-s±i5±i-r^ N _6-i >N _s-i ^^Jii "jv-i^-.JV-^i 

Finally the last four terms in brackets may be simplified (with most terms cancelling) leaving 27^, , . Again we 

arrive at 

/V 2 

b - 4 J ^' 
proving the conjecture. 
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